Performance of Block Jacobi-Davidson Eigensolvers
نویسندگان
چکیده
Jacobi-Davidson methods can efficiently compute a few eigenpairs of a large sparse matrix. Block variants of JacobiDavidson are known to be more robust than the standard algorithm, but they are usually avoided as the total number of floating point operations increases. We present the implementation of a block Jacobi-Davidson solver and show by detailed performance engineering and numerical experiments that the increase in operations is typically more than compensated by performance gains on modern architectures, giving a method that is both more efficient and robust than its single vector counterpart.
منابع مشابه
Performance evaluation of eigensolvers in nanostructure computations
We are concerned with the computation of electronic and optical properties of quantum dots. Using the Energy SCAN (ESCAN) method with empirical pseudopotentials, we compute interior eigenstates around the band gap which determine their properties. Numerically, this interior Hermitian eigenvalue problem poses several challenges, both with respect to accuracy and efficiency. Using these criteria,...
متن کاملInterior state computation of nano structures
We are concerned with the computation of electronic and optical properties of quantum dots. Using the Energy SCAN (ESCAN) method with empirical pseudopotentials, we compute interior eigenstates around the band gap which determine their properties. Numerically, this interior Hermitian eigenvalue problem poses several challenges, both with respect to accuracy and efficiency. Using these criteria,...
متن کاملThe Jacobi–Davidson method
The Jacobi–Davidson method is a popular technique to compute a few eigenpairs of large sparse matrices. Its introduction, about a decade ago, was motivated by the fact that standard eigensolvers often require an expensive factorization of the matrix to compute interior eigenvalues. Such a factorization may be infeasible for large matrices as arise in today’s large-scale simulations. In the Jaco...
متن کاملA Parallel Implementation of the Jacobi-Davidson Eigensolver and Its Application in a Plasma Turbulence Code
In the numerical solution of large-scale eigenvalue problems, Davidson-type methods are an increasingly popular alternative to Krylov eigensolvers. The main motivation is to avoid the expensive factorizations that are often needed by Krylov solvers when the problem is generalized or interior eigenvalues are desired. In Davidson-type methods, the factorization is replaced by iterative linear sol...
متن کاملStudies on Jacobi-Davidson, Rayleigh quotient iteration, inverse iteration generalized Davidson and Newton updates
We study Davidson-type subspace eigensolvers. Correction equations of Jacobi-Davidson and several other schemes are reviewed. New correction equations are derived. A general correction equation is constructed, existing correction equations may be considered as special cases of this general equation. The main theme of this study is to identify the essential common ingredient that leads to the ef...
متن کامل